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Abstract 

Fluctuations of the number of particles for the dilute interacting gas with Bose- 
Einstein condensate are considered. It is shown that in the Bogolubov theory these 
fluctuations are normal. The fluctuations of condensed as well as noncondensed 
particles are also normal both in canonical and grand canonical ensembles. 
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1 Introduction 



Fluctuations in systems with Bose-Einstein condensate have long been a topic of special 
attention, as can be inferred from reviews [1-3]. This problem plays an important role 
in the general context of quantum statistical mechanics. Interest in this problem has 
quickened in the past years due to intensive experimental and theoretical investigations 
of Bose-Einstein condensed atomic gases (see reviews [4-8]). 

The number-of-particle fluctuations in a uniform ideal Bose gas with Bose-Einstein 
condensate are known [1,2] to be anomalous, having the dispersion A 2 (N) ~ iV 2 , with 
the power of N larger than one. This anomalous behaviour comes from the fluctuations 
of the condensate, the corresponding dispersion being A 2 (A r o) ~ N 2 . The fluctuations of 
the noncondensed particles are also anomalous, with the dispersion A 2 (iVi) ~ iV 4 / 3 . Such 
anomalous fluctuations manifest the instability of the Bose-condensed ideal gas [3]. 

The problem of fluctuations in interacting Bose gases has been the arena of numerous 
discussions, with controversial conclusions (see review [3]). In the present paper, this 
problem is analysed in the frame of the Bogolubov theory [9-1 1] and it is shown that no 
anomalous fluctuations arise, but all fluctuations are normal. 



2 Number-of-particle fluctuations 

The measure of fluctuations for the number-of-particle operator N is provided by the 
dispersion 

A 2 (iV) = < N 2 > - < N > 2 , (1) 

where the angle brackets mean, as usual, statistical averaging. Dispersion (1) can be 
considered as an observable quantity, since it is directly connected with the isothermal 
compressibility, sound velocity, and structural factor, which are measurable quantities. In 
order to stress that the linkage between dispersion (1) and the observable quantities is 
general and exact, it is worth demonstrating how these relations can be rigorously derived. 

For an equilibrium system of N particles in volume V, with temperature T, it is easy 
to check that 

where ks is the Boltzmann constant, N =< iV >, and fi is chemical potential. The Gibbs 
potential G = G(T, P, N) is a function of temperature T, pressure P, and the number of 
particles N. From the differential 

dG = -S dT + V dP + n dN , 

taking into account that G = fiN, one has 

N dfi = —S dT + V dP , 

where S is entropy. The latter equality gives 

dfi \ 1 fdP\ 
dN) TV ~~p\dN) TV ' 
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where p = N/V is density, or, equivalently, 

(t) = 

\dfiJ TV \dP) TV 
For the variables N, P, and V, because of the differential 



there exists the equality 



'dN\ idP_\ fdV\ 
dP)v{dv) N \dNJp 



Using the Maxwell relation 

(dV\ fdp\ 
{dN) Tp \dP) TN 

and the property 

(d/A 1_ (dG\ 1 
\dP) TN -N\dP) TN -~p 

one gets 

(dN\ _ (dV\ 
[dp) TV ~~ P [dp) TN ' 
The definition of the isothermal compressibility 



owing to the above relations, can be rewritten as 



1 (<M\ J_ (dN\ , v 

Comparing Eqs. (4) with (2) yields 

(5) 

T Npk B T K ' 

In turn, the isothermal sound velocity s can be expressed through the compressibility as 

- 1 r.'.'i w 



m 



TN m P K T 



where m is particle mass. 

On the other hand, representing the operator N = J ^(r)tp(r) dr through the field 
operators ip\v) and i/)(r), from definition (1) it follows 

A 2 (N) = N + J p(r)p(r')[g(r, r) - I] drdr' , (7) 



which should be compared with the central structure factor 



S(Q) = 1 + 1 J p(r)p(r')[<?(r,r') - 1] dr dr' , (8) 

where p(r) =< ^'(r)^(r) >, and the pair correlation function is 

C V>t(r)V>t(OV>(QV>(r) > 
p(r)p(r') 

In this way, there exist exact relations 

A 2 (A) = Nk B T P K T = A = NS(0) , (10) 

valid for any equilibrium system, whether it is uniform or not. A stable equilibrium 
system requires that its compressibility be positive and finite, < kt < oo. An infinite 
compressibility would mean that the system immediately collapses or explodes. In a stable 
system, according to relations (10), the sound velocity and the central structure factor are 
also finite and positive. Thus, relations (10) tell us that in a stable equilibrium system 
the dispersion A 2 (A) must satisfy the stability condition 

A 2 (N) 

< ^ ' < OO (0 < K T < oo) (11) 

for any A > 0, including the thermodynamic limit N — > oo. Fluctuations satisfying the 
stability condition (11) are termed normal, while if condition (11) is not valid, fluctuations 
are called anomalous. Clearly, anomalous fluctuations imply instability. For instance, the 
ideal uniform Bose gas is unstable, since it has anomalous fluctuations with A 2 (A) ~ A 2 
and a divergent compressibility k t ~ A. Hydrodynamic equations for this gas with 
condensate are plagued by the appearance of unbound solutions [12]. 

Now let us pass to systems experiencing Bose-Einstein condensation. At the present 
time, a variety of trapped atomic gases is known to demonstrate this phenomenon (see 
recent review [3]). Several types of molecules, produced by means of Feshbach resonances, 
have been condensed. Bose-Einstein condensation might also appear in a system of semi- 
conductor biexcitons, consisting of two electrons plus two holes [13,14]. 

In the presence of condensate, the operator of the total number of particles is written 
as a sum 

A = A + A x , (12) 

whose terms correspond to condensed (A ) and noncondensed [N\] particles. Then dis- 
persion (1) takes the form 

A 2 (A) = A 2 (A ) + A 2 (Ai) + 2 cov(A , A\) , (13) 

in which the last term contains the covariance 

cov(A , N x ) = - < NqN x + AxAq > - < A >< N\ > . (14) 
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Since N and Ni are usually defined as commuting operators, covariance (14) reduces to 

cov(iVo,iVi) = < N Ni > - < N xN 1 > . (15) 

By the stability condition (11), we know that A 2 (iV) ~ N. Then the fundamental 
question is: Could the fluctuations of either condensed or noncondensed particles, or 
both, be anomalous, at the same time preserving the validity of the exact Eq. (13)? 
Recently there have appeared a number of papers stating that such fluctuations could 
be anomalous (see discussion in review [3]). In the following section we consider an 
interacting homogeneous system with Bose-Einstein condensate at low temperature and 
density, when the Bogolubov theory is applicable, and show that all fluctuations are 
normal. 



3 Dilute gas 

Current experiments with Bose-condensed atomic gases are typified by a rather low den- 
sity, such that pa^ <C 1, where a s is a scattering length. Atomic gases can be cooled down 
to very low temperatures, when practically all atoms are condensed, so that \Nq — N\ <C N. 
Under these conditions, the Bogolubov theory [9-11] becomes applicable. 

Let us start with the standard Hamiltonian for a uniform system of spinless atoms, 

H = J ^ f (r) ^- ^- V 2 - f?j i/>(r) dr+ \ J ^ t (r)^ t (r , )$(r - rXr')^(r) drdr' , (16) 

with an interaction potential $(r) = $(— r). For dilute gas, the latter is usually modelled 
by the contact interaction $(r) = (4irh 2 a s /m)5(r). However, for the sake of generality, 
we keep the form $(r). The chenical potential fi is included here in order to compare 
the results, corresponding to the grand canonical ensemble, with the original Bogolubov 
consideration [9-11] made in the frame of the canonical ensemble. 

Following the Bogolubov prescription let us separate the condensate by the shift in 
the field operator, 

^(r)=Vo + ^i(r), (17) 

where the condensate operator ipo does not depend on r owing to the system uniformity. 
The field operators of noncondensed particles, ipi(r), are assumed to possess the same 
Bose commutation relations as ip(r). From this it follows that the operator of condensed 
particles, ipo, commutes with all operators in the thermodynamic limit, when iV — > oo. 
The operators ip and ipi are assumed to be orthogonal, 

J ^i(r) dr = . 

Hence, N = N + Ni as in Eq. (12), with 

N = 4^j V , N, ee / Vl(r)Vi(r) dr . 
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For a uniform system, the field operators can be expanded in Fourier series 

k^O 



ipo = -7= , ^i(r) = a Wk{r) 



over the plane waves v?fc( r ) = e lk ' r /VV. Assuming that the interaction potential can also 
be Fourier expanded, we have 

$ ( r ) = \ E ^ r > = / $ (r)e- lk ' r dr . 

For the number-of-particle operators, we get 

N = 4«o , N x = a l a k ■ 

k^O 

Hamiltonian (16) is gauge invariant, hence the field operators are defined up to a global 
phase factor. Then one may chose a representation, where a is self-adjoint, such that 
a a = N and aja^ = N . Actually, this choice is absolutely not principal and is accepted 
with the sole aim to simplify the following formulas. 

After the Fourier transformation, Hamiltonian (16) acquires the form 

4 

H=Y,H {n \ (18) 

n=0 

in which the terms are grouped according to the number of noncondensed-particle oper- 
ators. In the zeroth order, 

where $0 = / 3>(r)dr. The first-order term vanishes, = 0. In the second order, 

(2 ( h 2 k 2 N \ t A> / + t t 

H = ~7f- + TT $ o - H I %a k + — 2^ §k [2ala k + a k aL k + a. k a k 

k^O \ / k^O 

where it is taken into account that = & k because of the symmetry of <&(r). The 
third-order term is 

^ (3) = 77 J2 ®i (aWk a k- q a q + a\a\_ q a k a<^ . 

V kq^Q 

Finally, the fourth-order term is written as 

H(4) = 7j\? 2 <$> q a\ala p+q a k - q . 

Remembering that we consider an almost condensed system, where Nq ~ N, the terms 
and # (4) can be treated as small perturbations to the terms up to the second order. 
Limiting ourselves by the latter terms, we have 
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+^T7 E ( 2a Iafc + aW-k + • (19) 

AV mo 

The contraction of the total fourth-order Hamiltonian (18) to its second-order part (19) 
is the first Bogolubov approximation. 

The next approximation is the replacement in Eq. (19) of the operator Nq by its 
average value Nq, which gives 

N 2 - - N 

H = w $0 ~ + S w *4 a k + S ( 2a fc° fc + 4 a -fc + a -fe a fe) > ( 20 ) 

where the notation 

^^ + f*°-" (21) 

is used. 

The Heisenberg equation of motion for ip is equivalent to the equation 

< —T- > = , 

which yields 

A* = $o + ^T E $ * ( 2 4a fc + 4«U + • (22) 

To remain in the frame of the second-order approximation in Hamiltonian (20), one has 
to use in Eq. (21) the zero-order part of the chemical potential (22), setting there 

Na 

v = Y % - (23) 

Another approximation is again based on the fact that iVo ~ N. Therefore one can 
replace in Eq. (20) iVo by N. This gives 

TV 2 N 
H = — $ + J2 Uka\a k — fiN + —— ^ ( 2 4«fc + a[al k + a- k a k ) . (24) 

Hamiltonian (24) is diagonalized by means of the Bogolubov canonical transformation 

a k = u k b k + v_ k b ] _ k , 

resulting in the diagonal form 

H = E + E e k blh - fiN , (25) 

in which the first term is the ground-state energy 

E = i Np% + i ]T (e k - u k - P § k ) , 

1 1 k^0 
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and the Bogolubov spectrum is 



e k = ^2p§ k u k + L) 2 k . 
The latter, with the notation for the sound velocity 



can be rewritten as 



e k = ^2mc 2 k u k + ul . 

Note that, due to Eqs. (21) and (23), the spectrum e k is gapless. The coefficient functions 
of the Bogolubov transformation are given by the equations 



yje 2 k + m?4 + e k 2 _ \je\ + m?c\ 



u k = 7T , % 



2e k ' 2e k 

In the case of the contact interaction, one has <3>fc = Aith 2 a s /m and c k = c= {h/m)^/Airpa s . 

The diagonal form (25) has been derived here starting with the Hamiltonian (16) 
corresponding to the grand canonical ensemble. The same Hamiltonian (25), up to the 
term —fiN, was derived by Bogolubov [9-11] in the frame of the canonical ensemble. 
Therefore all following calculations are actually identical for the grand canonical as well 
as canonical ensembles. 

A principal point must be stressed related to the transition from Eq. (19) to Eq. (20), 
when the operator N is replaced by its average value N . The equivalence of Hamiltonians 
(19) and (20) can be understood in two senses. In the strong sense, the equality of 
operators (19) and (20) requires that N = N . This, however, is not compulsory. And 
the equivalence of operators (19) and (20) can be understood in the weak sense, implying 
the equality of all their matrix elements. The latter can be reformulated as the equality on 
the average, so that the average values < H > for both forms (19) and (20) be coinciding. 
As is evident from expressions (19) and (20), their average values coincide then and only 
then, when the operator Nq is not correlated with the operators a\a k and a_ k a k . This 
means, in particular, the validity of the decoupling 

< iVoiVi > = < N >< Nx > . (26) 

Taking account of Eq. (26), we notice that covariance (15) vanishes, cov(A , A\) = 0. 
As a result, dispersion (13) becomes 

A 2 (iV) = A 2 (A> ) + A 2 (A\) . (27) 

All terms here are non-negative. Hence, if at least one of the dispersions, either A 2 (Ao) 
or A 2 (Ai), is anomalous, then the left-hand side, A 2 (A), is also anomalous. This, how- 
ever, would contradict the stability condition (11). Thus we come to an indispensable 
conclusion: In a stable equilibrium system, the fluctuations of condensed as well as of 
noncondensed atoms must be normal. 

If all fluctuations have to be normal, then how could one explain the appearance of 
anomalous fluctuations in a number of recent calculations (see review [3]) accomplished 
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on the basis of the same Bogolubov approximations? To answer this question, let us 
recollect how such anomalous fluctuations arise. The standard origin of their appearance 
in calculations is as follows. One considers the dispersion A 2 (Aq), containing the four- 
operator terms < a\a k a q a q > or < b k b k b q b q >. These are decoupled as < b k b k >< b q b q > 
+ < b\b q >< b k b\ >. Calculating A 2 (Aq), one meets the integral A 2 (Aq) ~ N / dk/k 2 . 
To avoid here the infrared divergence, one can replace the integral by a sum over the 
discretized spectrum of collective excitations [1,15]. This yields A 2 (Aq) ~ iV 4 / 3 . Another 
way could be by limiting the integration from below by k m i n m 1/L, with L ~ N 1 ^ 3 . 
Then again A 2 (Aq) ~ iV 4 / 3 in either canonical or grand canonical ensemble. This means, 
according to equality (27), that the dispersion A 2 (N) ~ iV 4 / 3 is anomalous, hence the 
system is unstable. 

In order to stress that the same type of the anomalous dispersion A 2 (Aq) arises in the 
grand canonical as well as in the canonical ensembles, let us consider A 2 (Aq) in the grand 
canonical ensemble, when it can be represented as 



A 2 (Aq) =k B T-^ <N X > - cov(A> , Aq) 

Since in the Bogolubov approximation cov(N , Aq) = 0, one has 

A 2 (iV 1 ) = k B T-^ <N,> . 



Keeping in mind that V- k = v k , the operator Aq can be written as 

= E [{ u t + v l) b lh + v\ + u k v k (b\bl k + bl k b k ) 

fc^O 



From here one gets 



< Ni > 
V 



-II 



2^4 



e% + m z c' 



coth 



2k B T 



dk 

J2^ 



For the case of the contact interaction, when c k = c, and involving the relation 



oj k = \j e\ + m 2 c 4 — mc 2 



we come to the form 

< Ni > V2 /mcV 
V ~ (2vr) 2 \~h) J XII 

in which 



+ x 2 -l-5 



1/2 



/ mc \ x 

coth — — — x , 

\2k B T I v^T^ 



dx , 



x = J5(2 + 5) , 5 = — 

v mc 2 

At zero temperature, one has the known result 

< Aq > 8 



N ~~ 3^ 



'pal (T = 0) . 



9 



For any nonzero temperatures, taking the derivative d < N\ > /dp and setting p = po&o, 
one recovers the same anomalous behaviour A 2 (A r 1 ) ~ iV / dx/x 2 ~ A 4 / 3 . 

However this anomaly is spurious and arises owing to a not self-consistent calculational 
procedure. Really, the Bogolubov approach is based on the Hamiltonian (19) of the 
second-order with respect to noncondensed-particle operators a k . All higher-order terms 
have been omitted from the initial Hamiltonian (18). Hence, such higher-order terms 
should also be omitted from the calculations of any other physical quantities. Calculating 
the fourth-order expression < A 2 > with the second-order approximation is not self- 
consistent, hence, is not correct. 

The total dispersion A 2 (A) can be found from Eq. (7). In order to be self-consistent 
in defining the pair correlation function g(r — r') = g(r, r'), one has to omit in the latter 
all terms of the order higher than two. Then one finds 



g(r) = 1 + - J (< a\a k > + < a k a_ k >) e ikr —- . (28) 



From here, 

/2 

[g(r) - 1] dr = - lim (< a[a k > + < a k a- k >) . 

p k^O \ ' 

This gives for the structural factor (8) 

5(0) = ^ (29) 
mc z V V m J 

And we find the dispersion (10) as 

A 2 (iV)=iV^4, (30) 

which is, of course, normal. If in the pair correlation function (28) we would retain the 
terms of the orders higher than the second with respect to a k , we would again get an 
anomalous dispersion. This, however, would not involve any senseful physics, but would 
simply mean the inconsistency of the calculations. 

The dispersion A 2 (A"i), being a fourth-order expression with respect to a k) is not a 
well-defined quantity in the frame of the second-order Bogolubov theory. For its defini- 
tion, it requires to invoke additional assumptions. Thus, we may follow the Bogolubov 
approximation, treating N — N as a classical quantity, because of which A 2 (A"o) = 0. 
Then in the grand canonical ensemble, 

A 2 (iV 1 ) = A 2 (iV) = N ^ (grand) . 

mc 2 

In the canonical ensemble, when N = N is fixed, we get 

A 2 (A> ) = A 2 (Ai) = N (canonical) , 

mc 1 

which results from the fact that A 2 (Ai) is the same in both ensembles. 
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Passing to the ideal gas, one should consider the limit $(r) — > 0, that is, c — > 0. Then 
A 2 (iV) diverges, as it should be for the ideal gas, which is unstable. The character of this 
divergence can be estimated as follows. The vanishing of $(r) can be accepted as being 
analogous to its diminishing at large distance, where it has the standard Lennard- Jones 
behaviour $(r) ~ 1/r 6 . In other words, $(r) ~ 1/iV 2 . Consequently, $ ~ ^/N, hence 
c ~ \/y/N. Then dispersion (30) diverges as A 2 (iV) ~ N 2 , which is typical of the ideal 
uniform Bose gas [1,2]. Contrary to the unstable Bose-condensed ideal gas, the interacting 
gas is stable, possessing always only normal fluctuations. 

The conclusion on the absence of anomalous fluctuations can be generalized for any 
stable equilibrium systems by rigorously proving the following theorem. If an observable 
quantity is represented as a sum A + B of two linearly independent terms, then the total 
dispersion 

A 2 (i + B) = A 2 (A) + A 2 (5) + 2 cov(A, B) , 

with the covariance 

cov(i, B) = \ <AB + BA> - <A><B>, 

is normal then and only then, when both partial dispersions A 2 (A) and A 2 (B) are normal. 
The proof will be presented in a separate publication. 

Acknowledgement. I am grateful for useful discussions to L.P. Pitaevskii and for 
helpful comments to Z. Idziaszek and E.P. Yukalova. The Mercator Professorship of the 
German Research Foundation is appreciated. 



11 



References 

[1] R.M. Ziff, G.E. Uhlenbeck, and M. Kac, Phys. Rep 32, 169 (1977). 

[2] D. ter Haar, Lectures on Selected Topics in Statistical Mechanics (Pergamon, Oxford, 
1977). 

[3] V.I. Yukalov, Laser Phys. Lett. 1, 435 (2004). 

[4] A.S. Parkins and D.F. Walls, Phys. Rep. 303, 1 (1998). 

[5] PW. Courteille, V.S. Bagnato, and VI. Yukalov, Laser Phys. 11, 659 (2001). 

[6] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Clarendon, Oxford, 
2003). 

[7] K. Bongs and K. Sengstock, Rep. Prog. Phys. 67, 907 (2004). 

[8] J.O. Andersen, Rev. Mod. Phys. 76, 599 (2004). 

[9] N.N. Bogolubov, J. Phys. (Moscow) 11, 23 (1947). 

[10] N.N. Bogolubov, Lectures on Quantum Statistics, Vol. 1 (Gordon and Breach, New- 
York, 1967). 

[11] N.N. Bogolubov, Lectures on Quantum Statistics, Vol. 2 (Gordon and Breach, New 
York, 1970). 

[12] V.I. Tsurkov, Majorant Catastrophe in Eulerian Gas Dynamic Equations for Bosons 
(Felicity, Whitefield, 1998). 

[13] D.W. Snoke, J.P Wolfe, and A. Mysyrowicz, Phys. Rev. B 41, 11171 (1990). 

[14] M. Hasuo, N. Nagasawa, T. Itoh, and A. Mysyrowicz, Phys. Rev. Lett. 70, 1303 
(1993). 

[15] S. Giorgini, L.P Pitaevskii, and S. Stringari, Phys. Rev. Lett. 80, 5040 (1998). 



12 



